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PREFACE 


This  report  describes  the  analysis  techniques  being  used  to  inter¬ 
pret  the  Wideband  Satellite  spaced-receiver  data  that  have  been  routinely 
recorded  at  the  Poker  Flat  and  Kwajalein  tracking  stations.  From  the 
spaced-receiver  data,  the  spatial  structure  of  the  diffraction  pattern  and 
its  apparent  drift  can  be  determined.  Through  careful  modeling,  this  in¬ 
formation  can  be  used  to  determine  the  anisotropy,  effective  height,  and 
true  drift  of  the  ionospheric  irregularities. 

Preliminary  results  from  the  Poker  Flat  station  confirm  the  sheet¬ 
like  irregularity  structure  that  we  hypothesized  to  explain  localized 
scintillation  enhancements  at  the  point  where  the  propagation  path  inter¬ 
cepts  the  local  L-shell.  The  auroral  irregularity  drifts  show  large  im¬ 
pulsive  variations  which  appear  to  be  associated  with  active  aurora. 

The  high  anisotropies  of  the  equatorial  data  cannot  be  resolved 
with  the  900-m  maximum  baseline.  As  far  as  irregularity  drifts  are  con¬ 
cerned,  the  apparent  east-west  drifts  are  largely  insensitive  to  axial 
ratio  increases  beyond  40:1.  Thus,  irregularity  heights  and  true  east- 
west  drifts  can  be  determined.  The  preliminary  data  show  evidence  of 
the  expected  eastward  drifts. 

The  analysis  is  done  in  three  parts.  First,  the  cross-correlation 
functions  for  all  possible  receiver  pairs  are  measured  and  processed  to 
determine  the  anisotropy  of  the  diffraction  pattern  as  a  function  of  time 
for  an  entire  pass.  The  measured  anisotropy  parameters  are  then  used  to 
determine  the  relative  pattern  drift.  Finally,  model  computations  are 
used  to  reproduce  the  measured  anisotropy  and  apparent  pattern  drift, 
thereby  giving  an  irregularity  height.  Once  the  height  is  known,  the 
relative  drifts  can  be  converted  to  true  irregularity  drifts  by  removing 
the  known  satellite  contribution. 
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Data  reduction  for  the  first  part  of  the  analysis  requires  much 
computation,  which  is  made  practical  by  using  an  array  processor.  A 
large  number  of  passes  have  now  been  processed  through  the  first  part  of 
the  analysis,  and  we  are  developing  procedures  to  automate  the  second  and 
third  parts  of  the  analysis.  This  will  allow  accumulation  of  a  data  base 
of  irregularity  heights  and  drifts,  which  should  enhance  our  understanding 
of  instability  mechanisms,  particularly  at  the  equator. 

The  model  computations  illustrate  the  importance  of  allowing  for 
the  geometric  dependence  of  the  scintillation.  This  has  important  rami¬ 
fications  for  translating  the  Wideband  satellite  results  to  other  geome¬ 
tries,  for  example,  geosynchronous  orbit,  as  well  as  for  predicting  scin¬ 
tillation  effects  in  general. 
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I  INTRODUCTION 


In  this  report  we  first  describe  a  general  correlation  function 
analysis  technique  that  allows  us  to  estimate  simultaneously  the  aniso¬ 
tropy  and  the  mean  drift  of  a  radio-wave  diffraction  pattern.  We  then 
develop  a  diffraction  theory  that  relates  the  anisotropy  and  mean  drift 
of  the  diffraction  pattern  to  the  anisotropy  and  mean  drift  of  the  iono¬ 
spheric  irregularities.  The  analysis  has  been  developed  to  interpret 
Wideband  Satellite  spaced- receiver  measurements. 

Spaced-receiver  measurements  have  been  used  to  decades  to  deduce 
ionospheric  drifts.  The  analysis  is  motivated  by  the  principle  that 
pairs  of  probes  sampling  a  rigid  but  moving  medium  will  produce  similar 
signals  displaced  in  time  by  an  amount  that  depends  on  the  probe  separa¬ 
tion  and  on  the  drift  component  along  the  axis  of  the  probes.  The  method 
of  similar  fades  (Mitra,  1949)  is  a  direct  application  of  this  principle. 

In  most  applications,  however,  one  must  allow  for  a  random  diffrac¬ 
tion  pattern  that  undergoes  some  rearrangement  as  it  drifts.  Variations 
of  the  correlation  method  developed  by  Briggs  et  al.  (1950)  have  proved 
to  be  superior  in  such  cases.  Indeed,  the  original  method  of  Briggs  et 
al.  gives  a  measure  of  the  random  drift  component  as  well  as  the  mean 
drift.  However,  because  the  random  component  does  not  admit  a  concise 
physical  interpretation,  it  is  not  used  much  now  (Wright  and  Pitteway, 
1978). 

Correlation  function  analysis  methods  nearly  always  assume  that 
surfaces  of  constant  correlation  are  ellipsoidal  and  therefore  give  rise 
to  elliptical  contours  in  any  plane.  Various  schemes  have  been  devised 
to  estimate  the  parameters  that  characterize  the  correlation  ellipses. 


* 

References  are  listed  at  the  end  of  the  report. 
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The  papers  by  Phillips  and  Spencer  (1955),  Yerg  (1955),  and  Fedor  (1967) 
are  representative. 

However,  these  schemes  require  restrictive  assumptions  about  the 
detailed  shape  of  the  spatial  and/or  temporal  correlation  functions. 

More  recently,  Armstrong  and  Coles  (1972)  devised  a  scheme  that  requires 
no  a  priori  assumptions  about  the  shape  of  the  autocorrelation  function 
or  its  time  variation.  The  Armstrong-Coles  method  uses  the  intersection 
points  of  the  various  correlation  functions  obtained  from  different  re¬ 
ceiver  pairs. 

If  the  anisotropy  is  known,  the  intersection  points  provide  suffi¬ 
cient  information  to  determine  the  mean  pattern  drift.  Thus,  by  assuming 
an  isotropic  diffraction  pattern,  Armstrong  and  Coles  (1972)  were  able  to 
measure  the  solar  wind  drift  speed  from  interplanetary  scintillation  ob¬ 
servations.  They  noted,  however,  that  if  the  delays  to  the  maxima  of  the 
correlation  functions  were  also  used,  the  anisotropy  of  the  diffraction 
pattern  could  be  determined  independently. 

It  is  interesting  that  the  delay  parameter  originally  used  to  deter¬ 
mine  the  pattern  drift  now  becomes  the  critical  parameter  for  resolving 
the  anisotropy  of  the  diffraction  pattern.  This  fact  can  be  ascertained 
by  recognizing  that  any  change  in  perspective  producing  a  high  degree  of 
correlation  along  a  particular  baseline  will  be  interpreted  as  a  large 
increase  in  the  pattern  drift  if  one  simply  divides  the  delay  parameter 
into  the  baseline  length.  This  effect  is  prominent  in  scintillation  data 
from  low-orbiting  satellites,  which  are  our  primary  interest. 

In  Section  III  we  describe  a  general  version  of  the  Armonstrong-Coles 
method  that  allows  us  to  simultaneously  estimate  the  anisotopy  and  the 
mean  pattern  drift  from  a  set  of  spaced-receiver  scintillation  observa¬ 
tions.  However,  because  of  the  complicated  interplay  between  the  actual 
three-dimensional  irregularity  anisotropy  and  the  propagation  geometry, 
in  Section  II  we  develop  the  diffraction  theory  that  relates  the  average 
structure  of  the  diffraction  pattern  to  a  general  model  of  the  average 
irregularity  structure. 
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This  report  primarily  is  intended  to  provide  a  detailed  and  self- 
contained  description  of  the  spaced-receiver  analysis  procedures  we  have 
developed  to  analyze  the  Wideband  Satellite  spaced-receiver  data.  [The 
Wideband  Satellite  Experiment  is  described  in  Fremouw  et  al.  (1978).] 

We  demonstrate  the  method,  using  data  from  both  auroral  and  equatorial 
stations,  in  Section  IV.  The  results  are  of  considerable  interest  in 
their  own  right,  and  they  will  be  developed  fully  in  other  reports. 

The  preliminary  results  from  the  auroral  zone  data  cannot  be  ex¬ 
plained  in  terms  of  ionospheric  irregularities  that  exhibit  a  rod- like 
(axially  symmetric)  anisotropy.  We  have  therefore  postulated  a  sheet¬ 
like  structure  in  which  the  irregularities  have  a  high  degree  of  coherence 
transverse  to  the  magnetic  field  along  the  direction  of  constant  magnetic 
field  strength,  i.e.,  along  the  L-shell.  Excellent  results  have  been  ob¬ 
tained  by  using  this  model. 

The  equatorial  data  are  consistent  with  the  expected  highly  elon¬ 
gated  rod-like  striations  oriented  along  the  magnetic  field  direction. 
Indeed,  the  axial  ratio  is  not  reliably  measured  with  the  900-m  north- 
south  baseline.  Hence,  the  drift  rates  cannot  be  directly  measured. 
Nonetheless,  the  east-west  component  of  the  apparent  pattern  drift  is 
insensitive  to  the  precise  value  of  the  axial  ratio,  a,  when  a  ^  40. 

Thus,  model  computations  can  still  be  used  to  deduce  the  true  irregularity 
drifts  and  heights. 

Both  the  equatorial  and  auroral  examples  clearly  reveal  the  impor¬ 
tance  of  allowing  for  the  effects  of  pattern  anisotropy  in  estimating 
irregularity  drifts. 
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II  THE  ANISOTROPY  MODEL 


In  this  section  we  develop  a  model  from  which  the  anisotropy  of  the 
diffraction  pattern  in  a  radio-wave  transmission  experiment  can  be  calcu¬ 
lated.  The  geometry  is  shown  in  Figure  1. 

The  reference  coordinate  system  is  located  within  the  region  of  the 
scattering  irregularities.  The  polar  angles  0  and  cp  determine  the  direc¬ 
tion  of  the  principal  ray  path.  If  we  ignore  refractive  effects,  the 
principal  ray  path  lies  along  the  line-of-sight .  The  vector  1?  is  in¬ 
clined  at  an  angle  in  the  xy  plane  and  lies  along  the  local  magnetic 
field  vector  at  P.  We  let  the  z  axis  point  downward,  so  that  i|;  is  the 
local  dip  angle.  A  system  of  receivers  is  configured  arbitrarily  about 
the  reference  point  R. 

It  is  clear  from  Figure  1  that  we  must  develop  a  fully  three- 
dimensional  model  to  characterize  the  diffracted  wavefield  in  the  neigh¬ 
borhood  of  R.  To  this  end,  we  first  consider  the  spatial  autocorrelation 
function  (ACF)  of  phase  in  the  linear  ray-optics  approximation.  That  is, 
we  calculate  phase  by  integrating  the  refractive  index  along  the  direction 
of  the  principal  ray  path.  It  is  shown  in  Rino  and  Fremouw  (1977)  that, 
by  using  this  linear  ray  optics  approximation,  the  spatial  ACF  of  phase, 

R  (Ap,  Az;  At),  has  the  form: 

60 

R60(Ap,  Az;  At)  =  <602>/ftf(Aps  -  vg  At)]  ,  (1) 

2 

where  (60  )  is  the  phase  variance  and  /p(y )  gives  the  shape  of  the  ACF.  The 

2 

fact  that  in  a  power-law  environment  (60)  depends  on  the  data  interval 
used  (Rino  and  Matthews,  1978)  does  not  affect  the  anisotropy  considera¬ 
tions. 

The  anisotropy  is  characterized  by  the  quadratic  function 
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SOURCE 


FIGURE  1  GEOMETRY  FOR  GENERALIZED  SPACE-RECEIVER 
ANALYSIS 


CApg2  -  BAps  Apg  +  AAoJ 
x _ x  y _ £ 

(AC  -  B2/4) 


It  is  convenient  to  write  Eq.  (2)  in  matrix  form  as 

f2(Apg)  =  ApgT  £  Aps  , 


where  £  is  the  2X2  matrix: 


The  anisotropy  coefficients  A,  B,  and  C,  depend  on  the  propagation  angles 
0  and  cp  as  well  as  on  the  axial  ratios  and  orientation  angles  that  charac¬ 
terize  the  anisotropy  of  the  ionospheric  irregularities.  Complete  defini¬ 
tions  of  A,  B,  and  C  are  given  in  Appendix  A. 

The  Az  dependence  of  R„  is  contained  in  Ap  ,  which  is  defined  in 
A  60  s 

terms  of  a  unit  vector  a  lying  along  the  xz-plane  projection  of  the 
t*  *T 

propagation  vector  k  as 

Ap  =  Ap  -  tan  0  a  Az  .  (4) 

S  K_ 

T 

A  similar  definition  relates  the  relative  drift  v  to  v  ,  v  ,  and  v  . 

s  x  y  z 

The  relative  drift  contains  the  true  irregularity  drift  and  the  velocity 
of  P  induced  by  the  source  motion.  We  note  that  v^  is  the  zonal  compo¬ 
nent,  and 


COS 

+  v^  sin  i]f 

(5a) 

sin  | 

-  V  COS  llf 

X 

J 

(5b) 

where  V||  and  v  are  the  parallel  and  perpendicular  meridional  components, 
respectively. 

The  two-dimensional  spectral  density  function  (SDF)  that  corresponds 
to  Eq.  (1)  has  the  form 


$^(K)  =  rg2\2L  sec2  0  ab (AN2)  Q  (Ak^  +  B^Ky  +  CKy) 


(6) 


where  r  is  the  classical  electron  radius,  \  is  the  wavelength,  L  is  the 
6  2 

layer  thickness,  (ANg)  is  the  electron  density  variance,  and  a  and  b  are 
axial  ratios  that  characterize  the  degree  of  anisotropy  of  the  irregulari¬ 
ties  along  and  transverse  to  the  magnetic  field  respectively.  The  func¬ 
tion  Q(q)  gives  the  shape  of  the  three-dimensional  irregularity  SDF,  and 
^iy)  and  Q(q)  are  related  by  the  expression 


/?(y) 


qJ  (qy)  Q(q)dq 


// 


qQ(q)dq 


(7) 
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Thus,  to  the  extent  that  the  anisotropy  of  the  measured  phase  and 
intensity  correlation  functions  is  characterized  by  Eq.  (2)  the  data  in¬ 
terpretation  is  comparatively  straightforward.  To  investigate  the  valid¬ 
ity  of  this  assumption,  therefore,  we  first  consider  the  intensity  ACF 
under  conditions  of  weak  scatter.  It  is  shown  in  Rino  and  Fremouw  (1977) 
that  the  intensity  SDF,  $6i(k)>  is  given  by  the  simple  relation 

$6l(K)  =  4§^(k)  sin2  [h(K)  ,  (8) 

where 

—*  2.  2  2. 
h(K)  =  K  +  tan  ©( K  cos  cp  +  K  sin  cp)  .  (9) 

x  y 


Thus,  the  intensity  ACF  has  the  general  form 


RgI(Ap,  Az;  z) 


4r2  \  L  sec2  6  ab(AN2) 


XJ^QC^  +  fiKxKy  +  CKy]  sin2  [h(K) 


X  cos  [k  •  A"p  ] 


dK 

(2n)‘ 


(10) 


Because  of  the  anisotropic  diffraction  term  in  Eq.  (10),  the  aniso¬ 
tropy  of  R._(Ap,  Az;  z)  is  strictly  characterized  by  Eq.  (2)  only  in 
PI  2 

special  situations.  As  an  example,  for  isotropic  irregularities  Ate 

2  x 

+  Bk  tc  +  Ck  =  h(K)  [see  Appendix  A],  and  Eq.  (2)  characterizes  the 
x  y  y 

anisotropy  of  R^.  If  we  change  variables  in  Eq.  (10)  so  that 


q  =  (k  cos  cp  +  K  sin  cp)  sec  6 
xx  y 

and 

q  =  (-K  sin  cp  +  K  cos  cp)  , 
y  x  y 


(11a) 


(lib) 


then 
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R6I(APS,;  Az>  z)  =  4re 


X^(L  sec  9)  ab(AN^) 


xff Q(A'q* +  B'q*qy  +  c/qy}  sin2  fq2  x'''z"'4f"~] 


X  cos  (q  •  Ap7)- 


_d^_ 


(2*)‘ 


(12) 


where 


Ap  7 

X 


Ap  7 


[cos  cp  cos  0  Ap  -  sin  9  Ap  ] 
s  s 

X  y 


[  sin  cp  cos  0  Ap  +  cos  cp  Ap  ] 
x  y 


(13a) 

(13b) 


The  coefficients  A7,  B7,  and  C7  are  defined  in  Rino  and  Fremouw  (1977). 

The  diffraction  term  in  Eq.  (12)  is  now  isotropic,  and  it  follows 

that  the  anisotropy  of  in  the  A^g  7  plane  is  given  by  Eq.  (2)  with  A, 

B,  and  C  replaced  by  A7,  B7,  and  C7.  Moreover,  the  same  transformation 

can  be  applied  to  the  phase,  so  that  both  the  phase  and  amplitude  have 

the  same  anisotropy  in  the  Apg7  plane  which  is  normal  to  the  propagation 
direction.  This  reference  system,  which  was  introduced  by  Briggs  et  al., 
(1950)  and  used  by  Singleton  (1970),  Moorcroft  and  Arima  (1972),  and 
others,  generally  simplifies  the  form  of  the  results,  but  it  is  not  com¬ 
pletely  general  unless  the  transrormation  defined  by  Eq.  (13)  is  used. 

The  analysis  presented  here  shows  that  there  are  subtle,  diffraction- 
induced  differences  between  the  anisotropy  of  the  spatial  structure  of 
diffracted  wavefield  and  that  given  by  Eq.  (2).  However,  for  moderate 
zenith  angles,  the  differences  are  quite  small.  If  we  consider  that 
multiple  scatter  effects  can  be  accommodated  by  iteratively  applying  the 
weak-scatter  results  (Rino,  1978),  it  follows  that  multiple-scatter  ef¬ 
fects  do  not  change  the  anisotropy,  although  the  shape  of  the  ACFs  of 
phase  and  intensity  may  change  dramatically  from  that  given  by  Eqs.  (1) 
and  (12)  respectively. 

One  must  also  keep  in  mind  that  in  a  real  (curved)  ionosphere  the 
angles  0,  cp,  and  change  along  the  propagation  path.  Thus,  for 


14 
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propagation  at  oblique  angles  in  an  extended  region,  there  is  a  geometric 
smearing  effect.  In  our  analysis,  we  ignore  the  geometric  smearing  and 
are  content  with  obtaining  a  set  of  effective  or  average  values  for  the 
anisotropy,  height,  and  drift  of  the  irregularities. 

Thus,  we  proceed  under  the  assumption  that  Eq.  (2)  adequately  charac¬ 
terizes  the  anisotropy  of  the  wavefield.  The  receiver  baselines  are 

specified  in  an  earth-fixed  geographic  coordinate  system  (Figure  2).  If 
-*T 

AR  =  (Ax  ,  Ay  ,  Az  )  represents  a  particular  baseline,  then  there  is  a 
K  K  K 

3X2  matrix  S'  (defined  in  Appendix  B)  such  that 

Aps  =  S'  .  •  (14) 

Hence, 

f2(Aps)  =  AlT  S/T£  S'  AR  .  (15) 

Ideally  one  would  determine  the  matrix  S/  and  then  use  Eq.  (14)  to 
calculate  Apg  to  interpret  a  particular  correlation  measurement.  As  we 
have  noted,  however,  /S/  depends  on  the  irregularity  height,  which  is  not 
known  a  priori.  As  we  shall  see  in  Section  III,  estimating  the  anisotropy 
and  drift  from  a  set  of  correlation  function  observations  requires  ex¬ 
tensive  computation.  Hence,  it  is  not  practical  to  iterate  the  entire 
procedure  for  a  number  of  heights. 

Thus,  we  have  adopted  an  approximate  procedure  by  ignoring  any 
height  differences  among  the  various  receivers.  We  take  £"  as  the  2X2 
matrix  obtained  by  deleting  the  s'^  and  s'^  elements  in  S/.  It  follows 
that 

A 

f2(Ap)  =  Tr'T  e  AR'  ,  (16) 

s  ^ 

where 

£  =  i/,T£.s"  (i7) 

— *  — » 

and  AR'  is  obtained  from  AR  by  deleting  the  Az  component. 

R 
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Similarly, 


C.E. 

FIGURE  2  RELATION  BETWEEN  LOCAL  GEOMAGNETIC  COORDINATE  SYSTEM  P 
AND  RECEIVER  SYSTEM  R 


The  anisotropy  and  drift  parameters  £,  and  vg  are  independent  of 
height.  The  fundamental  quantities  £  and  vg  can  then  be  computed  easily 
for  a  variety  of  heights  by  inverting  Eqs.  (17)  and  (18).  The  final 

A  A 

parameter  set  is  determined  by  calculating  the  variation  of  £  and  v  for 

~  s 

a  satellite  pass  with  an  assumed  set  of  axial  ratios  and  an  irregularity 
height.  The  theoretically  calculated  values  are  compared  to  the  measured 
values  to  obtain  a  best  fit. 


A  — f 

In  Section  III  we  shall  describe  a  procedure  for  estimating  £  and  vg 
from  a  set  of  spaced-receiver  correlation  function  measurements.  For 
display  purposes,  the  anisotropy  of  the  diffraction  pattern  is  charac¬ 
terized  by  the  axial  ratio, 


AR  = 


A  +  C  +  D 

AAA 

A  +  C  -  D 


1/2 


where 


A  AO 

D  =  [(A  -  C)Z 


*9  1/2 

+  B  ] 


(19a) 


(19b) 


and  by  the  orientation  angle  of  the  principal  irregularity  axis  relative 
to  the  y  axis. 


6 


1 

2 


atan 


(20) 
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Ill  SPACED- RECEIVER  ANALYSIS  FOR  ANISOTROPY  AND  DRIFT 


For  a  system  of  n  spaced  receivers,  there  are  n  temporal  correla¬ 
tion  functions,  R. .(At),  where  the  subscripts  i  and  j  denote  the  particu¬ 
lar  receivers  shown  in  Figure  1.  Based  on  the  model  described  in  Section 
II,  these  correlation  functions  will  have  the  general  form 


R.^At)  =  R  [f(2Pij  -  vgAt);  At]  .  (21) 

In  Eq.  (21)  we  allow  for  temporal  changes  in  the  diffraction  pattern  as 

it  drifts  with  the  apparent  velocity,  v  .  Thus,  R  .  need  not  equal  R... 

s  ii  j  j 

However,  for  i  ^  j,  R_  (At)  =  R^^-At),  so  only  n(n  +  1 )  /  2  correlation 
functions  need  actually  be  computed. 

2 

Let  us  now  assume  that  the  n  correlation  functions  are  numbered 
consecutively  from  1  to  n^  and  denote  a  particular  correlation  function 
and/or  baseline  in  this  set  by  a  single  superscript,  for  example 

R(j)(At)  =  R(j)[f(Ap(j)  -  vgAt) ;  At]  .  (22) 


Following  the  method  of  Armstrong  and  Coles  (1972),  we  identify  the 
intercept  time  delay,  t^,  such  that 

R<J)[f(Ap(J>  -  Vjk>;  TJk]  -  R(fc)[£(4p<k)  -  7sTjk);  Tjfc]  .  (23) 


2  2  1  3 

There  are  potentially  (n  )  -  n  -  2n(n  -  1)  =  n(n  -  3n  +  2)  such  inter- 

2 

sections.  The  subtracted  terms  account  for  the  intersections  of  the  n 
correlation  functions  with  themselves  and  the  2n(n  -  1)  intersections  be¬ 
tween  the  correlation  functions  and  their  mirror  images.  These  intersec¬ 
tions  are  independent  of  the  correlation  functions  themselves  and  there¬ 
fore  provide  no  useful  information. 
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By  using  the  matrix  notation  introduced  in  Section  II,  it  follows 
that  Eq.  (23)  holds  only  if 

(Ap(j)  -  vs  Tjkfe(Ap(j)  -  vsTjk)  =  (Ap(k)  -  vsT.j£(Ap(k)  -  7sT.k)  , 


vIg,(Ap(j)  -  Ap(k))  =  7^-[Ap(j)  e^p(j)  -  Ap(k)  CAp(k)]  .  (24) 


If  £  is  known,  Eq.  (24)  constitutes  an  overdetermined  set  of  equations 
that  can  be  solved  using  least  squares  techniques  for  v  .  This  scheme 
was  used  by  Armstrong  and  Coles  (1972)  to  determine  solar  wind  velocities 
from  spaced-receiver  interplanetary  scintillation  measurements. 

For  ionospheric  observations,  the  anisotropy  of  the  diffraction  pat¬ 
tern  is  unknown  and  must  be  determined.  To  determine  the  anisotropy,  we 
introduce  another  set  of  equations  by  considering  the  time,  lag  t.,  such 

that  R^^Ap^^;  t.)  is  maximized.  It  follows  from  Eqs.  (21)  or  (22) 

J  2  -*  ( i )  — 

that  t.  must  minimize  the  quadratic  form  f  (Ap  J  -  v  At).  By  direct 
J  s 

computation,  it  can  be  shown  that  Tj  must  satisfy  the  relation 

vT  e  Ap(j)  =  t.  V2  ,  (25) 

s  ~  H  j  eff 

where 

veff  =  f(V  *  (26) 

[Note  that  the  phase  ACF  given  by  Eq.  1  has  the  form  <60  >  %  (vef£  At).] 

Now,  if  we  substitute  Eq.  (25)  into  Eq.  (24),  we  obtain  the  system 
of  equations 

ipU)T[£/%£f]4’<J>  -  4p(k)I[e/'4f3  4p(k)  -  2Tjk<T.  -  Tk>  ,  (27) 


2  2 

which  can  be  solved  for  [6/v^^].  Since  ve££  is  simply  a  scale  factor, 
the  anisotropy  of  the  diffraction  patterns  is  completely  determined  by 
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Eq.  (27).  Indeed,  if  we  divide  Eq.  (24)  by  v^^,  it  is  clear  that  know¬ 
ing  [£/Vgff]  *-s  sufficient  to  determine  vg. 

We  can  write  Eq.  (27)  in  a  more  compact  form  by  defining  the  X  3 
matrix 


with  elements 


9  = 


D(1>  D(1)  D(1) 

xx  xy  yy 


(N  )  (N  )  (N  ) 

D  D  D 

-  xx  xy  yy  — ' 


D(i>  =  Ap(^  -  AP(k)‘ 

XX  rX  X 


D(i)  =  Ap(k)Ap(j)  -  Ap(j)Ap(k) 
xy  rx  y  x  y 


D(i)  -  4p<J)2  -  4p(k>; 

yy  y  y 


(28) 


(29a) 


(29b) 

(29c) 


3 

where  j  and  k  take  on  all  =  n(n  -  3n  +  2)  admissible  pairs  of  inter¬ 
section  values.  We  also  define  the  N^-vector, 


(N  )  (N  ) 

2t..  (t. 

Jk  j 


(nt)  T 

\  )] 


(30) 


and  the  3-vector, 


X  =  [C/ve2ff,  B/vJffi  A/v2ff]T 


eff’ 


eff- 


(31) 


Then  Eq.  (27)  can  be  written  as 

=  T 

which  has  the  least  squares  solution 


(32) 
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Similarly,  if  we  define  the  N_  X  2  matrix 


where 


)  -  B/2(Ap. 


B/2(Ap 


and  the  N  vector  g-  with  elements 


,  we  can  rewrite  Eq 


The  least-squares  solution  of  Eq.  (37)  is 


Equations  (33)  and  (38)  give  the  anisotropy  and  relative  pattern  velocity 
deduced  from  a  set  of  spaced-receiver  temporal  correlation  function  mea¬ 
surements. 


To  review  the  procedure,  from  a  set  of  n  spaced-receiver  measurements 
we  generate  n  temporal  ACFs  and  n(n  -  1 ) / 2  temporal  cross-correlation 


functions.  This  is  most  efficiently  done  by  using  fast  Fourier  trans¬ 
form  (FFT)  techniques.  We  then  reflect  the  cross  correlation  functions 

2 

to  generate  the  basic  set  of  n  correlation  functions. 


We  next  determine  T.,  the  time  lag  at  which  each  of  the  n  correla¬ 
tion  functions  achieves  its  maximum  value.  For  the  ACFs,  t.  is,  of 

J 

course,  identically  zero.  Moreover,  the  T ^  values  for  sets  of  comple¬ 
mentary  cross  correlation  functions  are  negatives  of  one  another. 


Finally,  we  determine  the  N  =  n(n  -  3n  +  2)  intersections.  Here 

again,  only  half  of  the  N^.  values  need  actually  be  computed,  since 

intersections  between  complementary  pairs  of  cross  correlation  functions 

are  negatives  of  one  another.  In  general,  however,  not  all  possible 

intersections  can  be  computed  accurately.  Two  situations  arise.  First, 

if  the  correlation  time  is  long  compared  to  the  data  interval,  then  small 

errors  in  the  correlation  function  estimate  can  lead  to  large  errors  in 

T.,  .  Second,  at  the  other  extreme  of  short  correlation  times,  the 
jk 

intersections  fall  into  the  noise. 


Thus,  in  ol  data  reduction  procedure  we  automatically  reject  T 
values  that  occur  beyond  some  preselected  decorrelation  time,  typically 
the  time  to  70%  decorrelation.  To  avoid  inaccuracies  due  to  long  corre¬ 
lation  times,  we  reject  intersections  that  lie  outside  a  preselected  sub¬ 
interval  within  the  data  interval.  One  can  increase  the  data  interval 
to  accommodate  long  correlation  times,  but  this  ultimately  compromises 
the  stationarity  of  the  data,  which  leads  to  poor  time  resolution. 

The  data  manipulations  we  have  described  here  are  rather  involved, 
but  they  are  practical  with  small,  general  purpose  computers  that  have 
array  processor  capability.  In  Section  IV,  we  present  examples  to  il¬ 
lustrate  the  technique  just  described. 
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IV  MODEL  COMPUTATIONS  AND  REPRESENTATIVE  EXAMPLES 


In  this  section  we  apply  the  theory  given  in  Section  II  and  the 
method  outlined  in  Section  III  to  spaced-receiver  data  from  the  wideband 
satellite  tracking  stations  at  Poker  Flat,  Alaska,  and  Kwajalein.  We 
first  consider  model  computations  based  on  the  theory,  use  them  to  verify 
the  data  analysis  method,  and  then  present  the  representative  Wideband 
data  sets.  Because  the  propagation  geometry  and  irregularity  structures 
at  high  latitudes  differ  so  from  those  at  the  equator,  the  data  from  each 
region  will  be  considered  separately. 

A.  High  Latitude  Data 

To  specify  the  propagation  angles  for  high  latitudes,  we  have  used 
a  representative  Wideband  Satellite  pass  over  the  Poker  Flat,  Alaska  re¬ 
ceiving  station.  The  satellite  is  in  a  1000-km  circular  sun-synchronous 
orbit  (Fremouw  et  al.,  1978).  The  nighttime  passes  progress  from  east 
(premidnight)  to  west  (postmidnight),  with  the  trajectory  of  the  iono¬ 
spheric  penetration  point  at  a  given  altitude  nearly  parallel  to  the 
local  geomagnetic  meridian. 

Figure  3  shows  the  ground  projection  of  the  350-km  penetration  point 
for  a  typical  postmidnight  pass.  The  ground  level  constant  magnetic 
L-shell  contours  are  also  shown  to  illustrate  the  near  meridional  align¬ 
ment  of  the  pass.  The  dip  angle  at  the  Poker  Flat  station  is  76°.  The 
propagation  vector  lies  within  the  L  shell  at  L  i  5.5.  For  an  overhead 
pass,  field-aligned  propagation  occurs  at  this  point. 

Now,  if  we  specify  the  axial  ratio  parameters  a,  b,  and  6  (see  Ap¬ 
pendix  A),  the  theory  developed  in  Section  II  can  be  used  to  compute  the 
axial  ratio  and  orientation  angle  of  the  diffraction  pattern  as  it  would 
be  measured  at  the  Poker  Flat  receiving  station.  Figure  4  shows  the  re¬ 
sults  of  such  computations  for  three  different  parameter  sets.  A 
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L  IONOSPHERIC 
1PENETRATION 
lpOINT  AT  350  km. 


POKER  FLAT* 


FIGURE  3  PATH  OF  350-km  PENETRATION  FOR  TYPICAL  POST-MIDNIGHT 
WIDEBAND  SATELLITE  PASS 


field-aligned  axially  symmetric  structure  (a  =  10,  b  =  1,  6  =  0)  produces 
highly  elongated  patterns  at  large  zenith  angles.  The  axial  ratio  de¬ 
creases  rapidly  toward  unity  as  the  propagation  vector  approaches  the 
magnetic  zenith.  The  orientation  of  the  principal  axis  measured  relative 
to  true  north  is  initially  along  the  geomagnetic  meridian  (S  s:  29°). 

Near  the  magnetic  zenith,  the  diffraction  pattern  undergoes  a  very  rapid 
change  in  orientation.  This  effect  is  easily  visualized  if  one  thinks 
of  the  shadow  of  a  rod  illuminated  by  a  moving  but  nonaligned  source. 

If  the  irregularities  have  a  high  degree  of  spatial  coherence  trans¬ 
verse  to  the  magnetic  meridian  plane  (b  >  1,  6  =  0)  as  well  as  along  the 
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FIGURE  4  AXIAL  RATIOS  AND  ORIENTATION  ANGLES  COMPUTED  FOR  THE 
HIGH-LATITUDE  WIDEBAND  PASS  SHOWN  IN  FIGURE  3,  USING 
DIFFERENT  IRREGULARITY  MODELS 


magnetic  field,  the  anisotropy  of  the  diffraction  pattern  changes  dra¬ 
matically.  As  can  be  seen  from  the  remaining  two  sets  of  curves  in 
Figure  4,  small  axial  ratios  now  occur  for  large  zenith  angles  while  the 
maximum  axial  ratio  is  achieved  near  the  field-aligned  direction  where 
the  diffraction  pattern  is  nearly  east-west  aligned.  Moreover,  the 
orientation  angle  variations  are  sensitive  to  the  degree  of  east-west 
anisotropy. 

We  see  from  these  examples  that,  at  high  latitudes,  the  anisotropy 
of  the  diffraction  pattern  is  a  very  sensitive  function  of  the  detailed 
irregularity  anisotropy.  To  show  this  in  another  way,  we  calculated  the 
delay  time  to  the  cross  correlation  peak  for  geomagnetic  east-west  and 
north-south  aligned  baselines. 
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The  actual  Wideband  Satellite  spaced-receiver  system  consists  of 
three  remote  antennas,  shown  schematically  in  Figure  5.  The  delay  param¬ 
eter  t  for  the  various  antenna  pairs,  is  easily  obtained  from  Eq.  (25). 

To  determine  v  ,  we  ignore  the  mean  irregularity  drift;  thus,  we  assume 
that  the  temporal  structure  of  the  diffraction  pattern  is  dictated  en¬ 
tirely  by  the  scan  velocity  induced  by  the  satellite  motion.  Since  the 
F-region  scan  velocity  is  not  less  than  ~  3  km/s,  this  is  a  reasonable 
assumption,  even  for  auroral  zone  phenomena. 

To  display  the  results,  we  convert  the  delay  parameter,  t  ,  to  an 
apparent  velocity  by  dividing  T\  into  the  corresponding  baseline  length. 
This  gives  the  traditional  apparent  velocity  parameter  as  defined  by 
Briggs  et  al.  (1950).  The  results  are  shown  in  Figure  6.  The  axially 
symmetric  model  produces  comparatively  small  and  fairly  uniform  apparent 
velocities  for  both  the  north-south  and  east-west  baselines. 

The  sheet  model,  by  comparison,  produces  an  extremely  large  apparent 
east-west  velocity  shear  at  the  point  where  the  propagation  path  lies 
within  the  local  L-shell.  Moreover,  the  apparent  north-south  velocity  is 
considerably  more  variable  than  the  corresponding  apparent  velocity  for 
axially  symmetric  irregularities.  These  results  illustrate  the  point 
made  in  Section  I;  the  delay  parameter  can  be  more  sensitive  to  the 
diffraction  pattern  anisotropy  than  to  the  relative  pattern  drift. 

For  comparison  with  the  apparent  pattern  velocities  shown  in  Figure  6, 

the  relative  pattern  drift,  vg,  for  this  same  pass  is  shown  in  Figure  7. 

In  this  case,  of  course,  the  relative  pattern  drift  is  due  totally  to 

satellite  motion.  The  velocities  shown  are  the  east-west  and  north-south 

components  of  v  .  One  should  keep  in  mind  that  v  is  influenced  by  the 
s  s 

vertical  drift  component  [see  Eqs.  (4)  and  (5)],  and  thus  there  will 
ultimately  be  some  ambiguity  in  determining  the  meridional  and  zonal  ir¬ 
regularity  drifts.  In  the  example,  the  drift  speed  is  minimum  at  ~  3  km/s 
at  the  point  of  closest  approach  and  increases  monotonically  for  xarger 
zenith  angles.  In  real  data,  departures  from  this  pattern  indicate  true 
irregularity  drifts  or  irregularity  height  changes. 
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FIGURE  5  WIDEBAND  SPACED-RECEIVER  GEOMETRY,  POKER  FLAT, 
ALASKA  AND  KWAJALEIN,  MARSHALL  ISLANDS 


Comparison  of  Figures  6  and  7  shows  clearly  that,  for  moving  beacons 
in  the  auroral  zone,  it  is  impossible  to  estimate  true  pattern  drifts, 
Particularly  east-west  drifts,  from  measured  apparent  velocities  without 
taking  the  irregularity  anisotropy  into  account. 

We  next  demonstrate  the  effectiveness  of  the  intersection  method, 
described  in  Section  III,  in  extracting  the  anisotropy  and  relative  drift 
from  the  actual  Wideband  system  under  ideal  conditions.  For  these 
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FIGURE  6  APPARENT  VELOCITIES  FOR  300-m  EAST-WEST  AND  900-m 

NORTH-SOUTH  BASELINES  FOR  THE  HIGH-LATITUDE  WIDEBAND 
PASS  SHOWN  IN  FIGURE  3  AND  THE  IRREGULARITY  MODELS 
USED  IN  FIGURE  4 


simulations,  the  time-changing  anisotropy  parameters  (axial  ratio,  orien¬ 
tation,  and  relative  pattern  velocities)  are  predicted  by  the  model--for 
each  antenna  and  between  each  of  the  antenna  pairs--and  then  used  in 
place  of  the  measured  correlation  functions  in  the  analysis  routine.  A 
simple  gaussian  form  for  R^\y)  [Eq.  (22)]  was  used,  with  a  decorrela¬ 
tion  time  similar  to  measured  values. 

The  satellite  pass  geometry  used  in  the  simulations  is  that  used  to 
produce  Figure  4.  For  the  three  anisotropy  cases  (sampled  every  20  s 
during  the  pass),  the  largest  discrepancies  between  the  input  and  output 
parameters  are:  axial  ratio  2%  and  orientation  0.01°.  The  errors  in 
extracting  the  relative  pattern  drifts  as  shown  in  Figure  7,  which  are 
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FIGURE  7  RELATIVE  IRREGULARITY  DRIFTS  FOR  THE  HIGH-LATITUDE 
WIDEBAND  PASS  SHOWN  IN  FIGURE  3,  SATELLITE-INDUCED 
COMPONENT  ONLY 


independent  of  the  irregularity  model  assumed,  are  better  than  1%,  except 
in  the  region  of  most  rapidly  changing  anisotropy,  where  errors  of  about 
10%  occur.  Overall  then,  under  ideal  conditions,  the  extraction  method 
works  very  well. 

u*  turn  now  to  actual  data.  Figure  8  shows  axial  ratio  and  orienta¬ 
tion  angle  estimates  obtained  by  applying  the  analysis  described  in 
Section  III  to  Wideband  Satellite  UHF  (413  MHz)  spaced-receiver  measure¬ 
ments  made  at  Poker  Flat  for  pass  number  12-11.  The  data  clearly  show 
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POKER  FLAT  12-11 
25  NOVEMBER  1976 


1130  1135  1140 

TIME,  UT 


FIGURE  8  MEASURED  AND  CALCULATED  AXIAL  RATIO  AND  ORIENTATION 
ANGLE.  POKER  FLAT  PASS  12-11,  25  NOVEMBER  1976 


the  enhanced  axial  ratios  near  the  center  of  the  pass  that  indicate  sheet¬ 
like  structures.  Indeed,  the  sheet-model  theory,  which  is  superimposed 
on  the  data  in  Figure  8,  fits  exceptionally  well  up  to  ~  1136  UT.  Beyond 
1136  UT,  the  UHF  perturbations  are  very  weak,  which  accounts  for  the  large 
data  scatter. 

As  one  should  expect,  the  model  anisotropy  computations  are  not  par¬ 
ticularly  height-sensitive.  On  the  other  hand,  the  magnitude  of  the  ap¬ 
parent  velocities  changes  significantly  with  height.  Figure  9  shows  the 
measured  apparent  velocities  for  the  same  data.  As  was  the  case  with 
Figure  6,  large  apparent  velocities  here  are  due  to  rapidly  changing 
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POKER  FLAT  12-11 
25  NOVEMBER  1976 


1130  1135  1140 


TIME,  UT 


FIGURE  9  MEASURED  AND  CALCULATED  APPARENT  VELOCITIES,  POKER 
FLAT  PASS  12-11,  25  NOVEMBER  1976 


irregularity  anisotropy,  with  a  prominent  shear  in  the  data  from  the 
east-west  baseline. 

The  sheet-model  data  are  shown  for  both  E-region  (110  km)  and  F- 
region  (350  km)  reference  altitudes.  The  measured  apparent  velocities 
generally  fall  between  the  two  model  calculations.  There  is,  however, 
a  localized  excursion  that  suggests  an  altitude  change  just  after  1132 
UT.  This  event  is  marked  by  the  arrows  in  Figure  9.  After  1133  UT,  the 
data  follow  the  E-region  curve. 

To  illustrate  the  computation  of  the  relative  velocity,  v  ,  we  used 
a  more  uniformly  disturbed  pass  than  12-11.  The  anisotropy  parameters 
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for  this  data  set,  pass  12-10  are  shown  in  Figure  10.  A  sheet  model  with 
b  <  a  seems  to  fit  the  data  better  than  the  model  with  a  =  b  =  10,  which 
was  used  for  pass  12-11.  The  corresponding  apparent  velocity  data  are 
shown  in  Figure  11.  Again,  the  apparent  velocity  shear  is  a  prominent 
feature  in  the  east-west  baseline  data. 

The  measured  relative  velocity  data  for  pass  12-10  are  shown  in 
Figure  12.  The  data  show  a  high  degree  of  variability,  which  may  be  real. 
However,  there  is  a  systematic  westward  enhancement  after  0947  UT,  which 
is  indicated  by  the  arrow  in  Figure  12.  Since  there  is  no  indication 
from  Figure  11  of  a  height  change  during  this  period,  the  excess  velocity 
can  be  attributed  to  irregularity  drifts.  Such  large  westward  velocities 
are  commonly  associated  with  premidnight  auroral  surges. 

Our  purpose  in  presenting  these  data  was  to  illustrate  the  technique 
of  simultaneously  estimating  the  anisotropy  and  relative  drifts  of  a 
radiowave  diffraction  pattern  in  an  extreme,  rapidly  changing  geometry. 

As  we  apply  the  analysis  to  the  large  body  of  acquired  data,  the  inter¬ 
pretation  will  be  developed  in  detail.  Nonetheless,  a  sheet-like  aniso¬ 
tropy  seems  to  be  generally  appropriate  for  auroral  zone  irregularities 
that  develop  within  the  region  of  the  continuous  (diffuse)  and  discrete 
aurora. 

B.  Equatorial  Data 

The  situation  at  the  equator  is  very  different  from  that  in  the 
auroral  zone,  both  in  terms  of  the  propagation  geometry  and  irregularity 
structure.  The  F-region  (350  km)  penetration  point  for  a  near-overhead 
pass  at  the  Wideband  site  on  Kwajalein,  Marshall  Islands  is  shown  in 
Figure  13.  Such  a  pass  scans  a  range  of  shallow  dip  angles,  from  ~  23° 
in  the  north  to  ~  -8°  in  the  south.  As  in  Alaska,  the  F-region  penetra¬ 
tion  scan  is  within  a  few  degrees  of  the  local  magnetic  meridian.  On 
the  ground,  the  spaced-receiver  arrangement  at  Kwajalein  is  as  shown  in 
Figure  5,  with  the  north-south  900-m  baseline  magnetically  aligned. 

Near  the  magnetic  equator,  the  high  conductivity  along  near-horizontal 
field  lines  creates  elongated  irregularity  structures.  Both  this 
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POKER  FLAT  12-10 
25  NOVEMBER  1976 


0945 


0950 


0955 


TIME,  UT 


FIGURE  10  MEASURED  AND  CALCULATED  AXIAL  RATIO  AND  ORIENTATION 
ANGLE,  POKER  FLAT  PASS  12-10,  25  NOVEMBER  1976 


elongation  and  the  near-field  alignment  of  the  satellite  scan  act  to  pro¬ 
duce  a  highly  elongated  diffraction  pattern  on  the  ground  for  the  entire 
satellite  pass.  To  show  this  quantitatively,  we  used  the  theory  given 
in  Section  II  to  compute  the  diffraction  pattern  parameters  for  a  near 
overhead  pass  at  Kwajalein.  The  ground  pattern  axial  ratio  and  orienta¬ 
tion  for  axially  symmetric  10:1  rods  (4  =  10,  b  =  1,  6  =  Q)  and  50:1  rods 
at  a  400-km  altitude  are  shown  in  Figure  14.  Projection  of  the  ionospheric 
structure  produces  a*ial  ratios  a  few  times  those  of  the  in-situ  irregu¬ 
larities  at  low  elevation  angles,  and  essentially  reproduces  the  in-situ 
anisitropy  near  the  closest  approach.  The  orientation  of  the  ground 
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FIGURE  11  MEASURED  AND  CALCULATED  APPARENT  VELOCITIES,  POKER 
FLAT  PASS  12-10,  25  NOVEMBER  1976 


correlation  ellipse  remains  near  10  to  12  deg  from  geographic  north 
during  the  entire  pass. 

The  apparent  velocities  for  this  same  pass  geometry  and  irregularity 
structure  are  shown  in  Figure  15.  The  east-west  apparent  velocity  for 
the  rod  model  is  a  few  hundred  meters  per  second  and  decreases  with  in¬ 
creasing  axial  ratio.  Beyond  an  axial  ratio  of  about  50:1,  the  east-west 
apparent  velocity  essentially  saturates  and  no  longer  decreases  with  in¬ 
creasing  axial  ratio.  Because  the  propagation  path  is  not  quite  in  the 
local  magnetic  meridian,  the  satellite  scan  produces  large,  north-south, 
apparent  velocity  shears.  The  increase  from  10:1  to  50:1  rods  dramatically 
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FIGURE  12  MEASURED  AND  CALCULATED  RELATIVE  VELOCITIES.  POKER 
FLAT  PASS  12-10,  25  NOVEMBER  1976 
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FIGURE  13  MAP  OF  THE  350-km  IONOSPHERIC  PENETRATION  POINT  FOR  A 
TYPICAL  HIGH-ELEVATION  WIDEBAND  PASS  AT  KWAJALEIN. 

Lines  of  constant  dip  latitude  are  superimposed 

affects  the  location  of  the  shear  early  in  the  pass,  and  actually  pro¬ 
duces  a  second  shear  near  the  end. 

For  comparison  with  these  apparent  velocity  drifts,  the  relative 
drift  due  to  satellite  scan  motion  is  shown  in  Figure  16.  Comparison  of 
these  relative  velocities  and  the  apparent  velocities  in  Figure  15  show 
again  (as  at  high  latitudes)  that  for  time-changing  propagation  geometry, 
removal  of  the  anisotropy  from  the  ground  diffraction  pattern  is  essential 
in  interpreting  apparent  pattern  velocities. 


36 


FIGURE  14  AXIAL  RATIOS  AND  ORIENTATION  ANGLES  COMPUTED  FOR 
THE  EQUATORIAL  WIDEBAND  PASS  SHOWN  IN  FIGURE  13 
AND  ROD-LIKE  IRREGULARITY  MODELS 


One  point  of  concern  in  using  the  intersection  method  on  the  equa¬ 
torial  data  is  whether  the  north-south  baseline  spacing  of  ~  900  m  is 
adequate  to  measure  the  extreme  pattern  anisotropy  that  may  occur.  To 
ensure  that  the  method  can,  in  fact,  resolve  large  axial  ratios,  a  series 
of  simulations,  similar  to  those  discussed  earlier,  have  been  carried 
out.  Gaussian  correlation  functions  with  decorrelation  times  comparable 
to  those  in  the  data  were  again  used  to  fully  simulate  the  time-changing 
ground  diffraction  pattern. 

Even  with  our  relatively  short  baseline,  the  intersection  method 
succeeds  in  extracting  the  correct  axial  ratio  for  the  rod  model  and 
ground-pattern  axial  ratios  up  to  about  200:1  with  small  errors.  For 
example,  100:1  in-situ  rods  produce  pattern  axial  ratios  up  to  about 
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FIGURE  15  APPARENT  VELOCITIES  FOR  300-m  EAST-WEST  AND  900-m  NORTH-SOUTH 
BASELINES  FOR  THE  WIDEBAND  EQUATORIAL  PASS  SHOWN  IN  FIGURE  13, 
AND  THE  IRREGULARITY  MODELS  USED  IN  FIGURE  14 

400:1  at  the  horizons.  Pattern  axial  ratios  up  to  about  200:1  are  ex¬ 
tracted  with  errors  less  than  1%;  between  200:1  and  400:1,  random  errors 
of  up  to  12%  are  observed.  It  is  important  to  note,  however,  that  al¬ 
though  errors  are  seen  for  these  high  axial  ratios,  the  method  never 

2 

fails  to  obtain  an  estimate  [i.e.,  AC  >  B  /4  for  all  cases,  see  Eq.  (3)]. 
In  no  case  is  the  orientation  error  greater  than  0,01°. 

The  analysis  method  does  not  do  as  well  in  extracting  the  relative 
drift  velocities.  For  pattern  axial  ratios  up  to  about  100:1,  the  errors 
in  the  extracted  east-west  and  north-south  relative  velocities  are  below 
~  3%.  For  axial  ratios  between  100:1  and  150:1,  random  errors  of  up  to 
about  20%  are  seen;  beyond  about  150:1,  the  intersection  method  generally 
fails  to  give  reasonable  results  for  relative  drift  velocities. 
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FIGURE  16  RELATIVE  IRREGULARITY  DRIFTS  FOR  THE  EQUATORIAL  WIDEBAND  PASS 
SHOWN  IN  FIGURE  14,  SATELLITE-INDUCED  COMPONENT  ONLY 


All  in  all,  this  indicates  that  under  ideal  conditions  we  should 
obtain  axial  ratio,  orientation,  and  relative  drifts  accurate  to  a  few 
percent  for  100:1  in  situ,  rod-like,  F-region  irregularities  for  more 
than  half  of  a  typical  satellite  pass.  For  the  remainder  of  the  pass 
good  anisotropy  but  poor  relative  drift  estimates  are  expected. 

We  next  consider  the  experimental  data  fcr  a  high-elevation  Kwajalein 
data  set,  pass  27-19.  The  geometry  of  this  pass  was  used  to  generate  the 
model  calculations  of  Figures  14  through  16.  The  north-south  and  east-west 
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apparent  velocities  for  pass  27-19  are  shown  in  Figure  17.  As  was  pointed 
out  in  conjunction  with  Figure  15,  the  axial  ratio  has  a  strong  effect  on 
the  location  of  apparent  north-south  velocity  shears.  In  this  case,  the 
location  of  the  shear  near  1219  UT  is  best  matched  by  in-situ  F-region 
rods  with  an  axial  ratio  of  ~  100:1,  the  signature  of  which  is  superim¬ 
posed  in  Figure  15.  Although  the  shear  location  is  well-matched  in  time 
in  this  comparison,  there  is  a  general  underestimation  of  both  the  north- 
south  and  east-west  apparent  velocities.  This  is  most  likely  an  indica¬ 
tion  of  true  irregularity  drift.  We  return  to  this  point  below. 

The  observed  axial  ratio  and  orientation  for  this  same  pass  are 
shown  in  Figure  18.  The  agreement  between  the  data  and  the  500-km  100:1 
rod  model  is  good  for  orientation  but  poor  for  axial  ratio.  The  good 
correspondence  between  the  model/data  orientation  is  not  surprising, 
since  the  correlation  ellipse  orientation  is  generally  insensitive  to 
changing  anisotropy  as  long  as  the  irregularities  remain  field-aligned. 

The  reason  that  the  measured  axial  ratio  is  so  low  relative  to  the 
model  calculations  is  more  difficult  to  explain.  First,  we  have  strong 
evidence  for  irregularity  axial  ratios  of  ~  100:1  in  the  apparent  velocity 
data  of  Figure  17.  Second,  we  have  shown  that  the  intersection  method 
extracts  axial  ratios  well,  under  ideal  conditions,  despite  the  extreme 
pattern  elongation.  Scrutiny  of  both  the  intersection  method  and  the 
data  indicates  that  one  reason  for  poor  performance  in  extracting  axial 
ratio  from  real  data  could  be  that  we  do  not  have  accurate  enough  esti¬ 
mates  for  the  spaced-receiver  antenna  locations.  This  and  other  possi¬ 
bilities  are  being  pursued  further. 

Although  we  currently  have  difficulty  in  obtaining  anisotropy  and 
relative  drifts  from  the  equatorial  data,  considerable  information  can  be 
obtained  from  the  apparent  velocity  data.  For  example,  Figure  19  compares 
the  measured  apparent  velocity  data  and  the  model  calculations  with  100  m/s 
200  m/s,  and  300  m/s  eastward  F-region  irregularity  drifts.  These  drift 
magnitudes  are  consistent  with  other  measures  of  the  nighttime  equatorial 
ionosphere  (e.g.,  those  of  Woodman  and  La  Hoz,  1976).  Both  the  north- 
south  and  east-west  measured  apparent  velocities  are  bracketed  well  by 


FIGURE  17  MEASURED  AND  CALCULATED  APPARENT  VELOCITIES,  KWAJALEIN  PASS 
27-19,  26  AUGUST  1977 


100:1  F-region  rods,  drifting  between  100  and  300  m/s  eastward.  The 
primary  shortcoming  of  this  approach  is  that  it  is  generally  insensitive 
to  any  north-south  drifts,  which  thus  cannot  be  extracted.  The  general 
ambiguity  between  vertical  and  horizontal  drifts  [see  Eqs.  (4)  and  (5)] 
still  remains. 

As  with  the  high-latitude  data,  a  large  number  of  equatorial  passes 
suitable  for  spaced-receiver  analysis  have  been  collected.  Detailed 
interpretation  of  this  data  and  continued  refinement  of  the  method  are 
under  way.  To  date,  information  such  as  that  in  Figures  14  through  19 
gives  us  preliminary  but  firm  indications  that  axially  symmetric,  highly 
elongated,  rod-like  irregularities  produce  the  bulk  of  the  equatorial 
scintillation. 
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AXIAL  RATIO  ORIENT A 


FIGURE  18  MEASURED  AND  CALCULATED  AXIAL  RATIO  AND  ORIENTATION  ANGLE, 
KWAJALEIN  PASS  27-19,  26  AUGUST  1977 
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V  DISCUSSION 


We  have  developed  and  demonstrated  a  technique  (actually  suggested 
by  Armstrong  and  Coles,  1972)  for  simultaneously  estimating  the  anisotropy 
and  relative  irregularity  drifts  of  a  radio  wave  diffraction  pattern. 
Because  of  the  complicated  geometric  dependence  of  the  scattering  phe¬ 
nomenon,  the  apparent  velocities  can  differ  considerably  from  Che  rela¬ 
tive  irregularity  drifts — thus  the  need  for  the  rather  elaborate  proce¬ 
dure. 

By  using  model  computations,  we  have  inferred  actual  irregularity 
anisotropy  parameters  and  heights  that  could  explain  the  data.  Moreover, 
by  calculating  the  satellite-induced  component  of  the  relative  irregu¬ 
larity  drifts  we  were  able  to  estimate  true  irregularity  drifts.  The 
preliminary  results  validate  the  technique  and  are  interesting  in  their 
own  right. 

Nonetheless,  one  would  like  some  guidelines  to  evaluate  the  uncer¬ 
tainty  in  the  inferred  parameters.  We  have  not  yet  attempted  to  derive 
error  bars.  In  the  analysis,  however,  we  used  the  magnitude  of  the  de¬ 
terminants  of  the  square  matrices  in  Eqs.  (33)  and  (38)  to  establish  an 
acceptability  threshold  for  individual  parameter  estimates.  Beyond  this, 
we  relied  on  the  agreement  between  the  model  calculations  and  the  param¬ 
eter  estimates  to  establish  confidence  in  the  results. 

There  are  also  practical  considerations.  The  sample  interval  and 
data  block  length  critically  influence  the  precision  that  can  be  obtained. 
The  Wideband  Satellite  spaced-receiver  system  uses  four  fully  coherent 
spaced  receivers.  Thus,  the  anisotropy  and  drift  parameters  can  be  inde¬ 
pendently  determined  from  the  phase  and  amplitude  correlation  functions. 
Since  phase  decorrelates  much  more  slowly  than  amplitude,  the  latter  is 
generally  easier  to  use. 
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The  amplitude  and  phase  data  are  detrended  as  described  in  Fremouw 
et  al.  (1978)  prior  to  the  correlation  function  analysis.  This  procedure 
effectively  establishes  a  common  reference  level  for  the  signals,  thereby 
minimizing  normalization  errors.  The.cross  correlation  delays  t.  and  the 
intersections  are  precisely  determined  by  polynomial  interpolation. 

By  using  an  array  processor,  the  data  reduction  for  an  entire  ~  15-min 
pass  is  performed  in  less  than  1  hour. 

In  summary,  the  apparent  velocity  parameter  cannot  be  simply  inter¬ 
preted  to  obtain  pattern  drifts.  Rather,  an  analysis  procedure  must  be 
employed  that  simultaneously  determines  the  anisotropy  and  relative  drift 
of  the  diffraction  pattern.  We  have  described  and  demonstrated  such  a 
procedure  in  this  report.  In  the  next  phase  of  analysis,  procedures  will 
be  developed  to  automate  the  extraction  of  height  and  true  drift  estimates 
from  the  measured  anisotropy  and  apparent  pattern  drifts. 
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Appendix  A 


ANISOTROPY  COEFFICIENTS 

In  Rino  and  Fremouw  (1977)  it  is  shown  that  the  two-dimensional 

phase  SDF,  $  Oic),  is  related  to  the  three-dimensional  irregularity  SDF 
60 

$AN  (K»KZ)  by  the  relation 
e 

*  0?)  =  r2\2L  sec2  0  $  (K,  -tan  0  a  K)  .  (A-l) 

60  e  ANg  kt 


Now,  if  we  assume  that  there  exists  a  coordinate  system  (r,  s,  t)  in 

which  the  three-dimensional  spatial  correlation  depends  only  on 

2  2  2  1/2 

[Ar  +  (As/a)  +  (At/r)  ]  ,  with  a  and  b  s  1,  the  corresponding 

functional  dependence  of  (K,  Kz)  on  K  and  Kz  can  be  computed  by 

0 

straightforward  manipulations. 

For  the  coordinate  system  of  Figure  1,  the  functional  dependence  is 

-*  X  A  -*  1/2  A 

given  in  terms  of  the  quadratic  form  [(K»K)  C(K,  K  )]  where  C  is 
a  3  X  3  matrix  with  elements: 


A  2  2  ,  .  2  .,2  .2  2  .  v 

C^  =  a  cos  |  +  sin  ^(b  sin  6  +  cos  6) 


a  2  2  2 

C  =  b  cos  6  +  sin  6 


A  2  2  2  2  2  2 
C^^  =  a  sin  ^  +  cos  ^(b  sin  6  +  cos  6) 


A  A  2 

C12  =  C21  =  (b  "  1)  sin  i}'  sin  6  cos  6 


A  A 


2  2  2  2 
C13  =  C31  =  “  b  sin  6  -  cos  6)  sin  ^  cos  i|i 


A  A 


C23  =  C32  =  "^b  ”  ^  cos  t  sin  ®  cos  6 


(A- 2a) 
(A- 2b) 
(A- 2c) 
(A- 2d) 
(A-2e) 
(A-2f) 
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Appendix  B 


SOME  USEFUL  COORDINATE  TRANSFORMATIONS 


Earth-fixed  spaced  receivers  are  most  conveniently  measured  in  a 

system  as  shown  in  Figure  2  where  z  is  positive  upward  and  y  is 

K  R 

positive  northward.  The  rotation  to  any  other  geographic  system,  say 
x  y  z  ,  is  given  by  the  3x3  matrix  R  with  elements: 


PP  P  ' 

ril  =  C°S  XP  C°S  XR  +  Sin  XP  Sin  ^R  (B-la) 
rl2  =  ^"cos  Xp  sin  XR  +  sin  XP  cos  XR)  sin  'Pr  (B-lb) 
r13  =  (cos  Xp  sin  XR  -  sin  Xp  cos  \R)  cos  cpR  (B-lc) 
r 21  ~  (-  sin  Xp  cos  \R  +  cos  Xp  sin  XR)  sin  cpp  (B-ld) 


*  (sin  Xp  sin  XR  +  cos  Xp  cos  XR)  sin  cpp  sin  cpR  +  cos  cpp  cos  cpR  (B-le) 

*  -(sin  Xp  sin  XR  +  cos  xp  cos  XR)  sin  cpp  cos  cpR  +  cos  cpp  sin  cpR  (B-lf) 
=  (sin  xp  cos  XR  -  sin  XR  cos  Xp)  cos  cpp  (B-lg) 
=  -(sin  Xp  sin  XR  +  cos  xp  cos  XR)  cos  cpp  sin  cpR  +  sin  cpp  cos  cpR  (B-lh) 
=  (sin  xp  sin  XR  +  cos  Xp  cos  XR)  cos  cpp  cos  cpR  +  sin  cpp  sin  cpR  (B-li) 


where  XR  and  xp  denote  the  respective  longitudes  and  cpR  denotes  the 
respective  latitudes. 

To  transform  a  point  from  the  x  y  z  system  to  the  xyz  system,  a 

K  K  K 

second  rotation  and  a  reflection  are  required;  this  is  implemented  by 
multiplying  by  the  matrix  D,  If  5  represents  the  magnetic  declination 
angle,  then 


51 


The  angle  i|r  is  shown  in  Figure  1.  The  angle  6  is  the  orientation  angle 
of  the  transverse  irregularity  axis.  Thus,  if  we  take  s  to  lie  along 
the  local  magnetic  field  direction,  \|/  is  the  dip  angle. 


The  coefficients  A,  B,  and  C  are  readily  deduced  by  substituting 

-tan  0  ak  •  K  for  K  and  collecting  the  terms  multiplying  k^,  K  tC  »  and 
2  T  2  x  x  y 

Ky.  The  result  of  these  manipulations  gives  the  desired  expressions: 


a  a  2  2  A 

A  =  [Cu  +  C33  tan  9  cos  cp  -  2C1  ^  tan  9  cos  cp] 


13 


(A-3a) 


B  =  2[C^2  +  ban  9  sin  cp  cos  cp 


-tan  9  (Cj^  cp  +  C?3  cos  cp)  ] 


23 


(A-3b) 


,  2 


C  =  [C22  +  C33  tan  9  sin  cp  -  2C23  tan  9  sin  cp]  .  (A-3c) 


We  note  that  if  a  =  b  =  1  (isotropic  irregularities),  then 


2  2 

A  =  1  +  tan  9  cos  cp 
2 

B  =  2  tan  9  sin  cp  cos  cp 
2  2 

C  =  1  +  tan  9  sin  cp 


(A-4a) 
(A-4b) 
(A-4c ) 


Thus,  the  diffraction  pattern  is  not,  in  general,  isotropic,  even  if  the 
irregularities  themselves  are  isotropic. 
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SOME  USEFUL  COORDINATE  TRANSFORMATIONS 

Earth-fixed  spaced  receivers  are  most  conveniently  measured  in  a 

system  as  shown  in  Figure  2  where  z  is  positive  upward  and  y  is 

K  R 

positive  northward.  The  rotation  to  any  other  geographic  system,  say 


x  y  z  ,  is  given  by  the  3x3  matrix  R  with  elements: 

P  P  P 

r i i  =  cos  Xp  cos  XR  +  sin  \p  sin  \R  (B-Ia) 

ri2  =  ("cos  xp  sin  XR  +  sin  Xp  cos  XR)  sin  'pR  (B-lb) 

rl3  =  (cos  Xp  Sin  XR  -  sin  Xp  cos  XR)  cos  cpR  (B-lc) 

r21  =  Sin  COS  XR  +  C0S  ^P  Sin  Sin  °Pp  (B-ld) 

r 22  =  (sin  Xp  sin  XR  +  cos  Xp  cos  XR)  sin  cpp  sin  cpR  +  cos  cpp  cos  cpR  (B-le) 


r23  =  “(sin  Xp  sin  XR  +  cos  Xp  cos  XR)  sin  cpp  cos  cpR  +  cos  9p  sin  cpR  (B-lf) 

r31  =  (sin  xp  cos  XR  -  sin  XR  cos  Xp)  cos  cpp  (B-lg) 

r32  =  -(sin  Xp  sin  XR  +  cos  xp  cos  XR)  cos  cpp  sin  cpR  +  sin  cpp  cos  cpR  (B-lh) 

r33  =  ^sin  XP  sin  XR  +  cos  ^p  cos  Xr)  cos  9p  cos  cpR  +  sin  cpp  sin  cpR  (B-li) 

where  XR  and  Xp  denote  the  respective  longitudes  and  ^>R  denotes  the 
respective  latitudes. 

To  transform  a  point  from  the  xRyRzR  system  to  the  xyz  system,  a 
second  rotation  and  a  reflection  are  required;  this  is  implemented  by 
multiplying  by  the  matrix  D.  If  g  represents  the  magnetic  declination 
angle,  then 
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(sin  6 
cos  § 
0 


cos  6 
-sin  6 
0 


(B-2) 


The  overall  transformation  is  denoted  by  the  3x3  matrix  S  =  DR  with 

/V  A/W 

elements  s .  . . 

ij 

Now,  as  shown  in  Section  III,  contours  of  constant  correlations 
of  intensity  or  phase  are  defined  by  constant  values  of 


f2(Aps)  =  ApsT  e  Aps 


(B-3) 


If  AR  represents  a  vector  in  the  xDyDz  system,  the  2-vector  Ap  is  ob- 

K  K  K  S 

tained  by  first  generating  £  Z\R  and  then  using  the  definition  of  Apg 
in  Eq.  (4).  The  overall  manipulation  can  be  represented  by  the  2X3 
matrix  S'  with  elements: 


s 


s 


s 


s 


s 


s 


11 

/ 

12 

/ 

13 

/ 

21 

/ 

22 

/ 

23 


s^  +  tan  0  cos  cp  r^ 
s12  +  tan  0  cos  cp  r^ 
s^^  +  tan  6  >;os  cp  r^^ 
s^  +  tan  9  cos  cp  r^ 
s22  +  tan  9  sin  cp  r32 

S23  +  tan  ®  sin  'P  r33 


In  terms  of  S',  therefore, 

f2(Aps)  =  2rT  S'T  £  £  AR 


(B-4a) 

(B-4b) 


(B-4c) 
(B-4d ) 
(B-4e) 
(B-4f ) 


(B-5) 
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Electrospace  Systems,  Inc. 

Johns  Hopkins  University 

ATTN 

H.  Logston 

Applied  Physics  Lab. 

ATTN 

Document  Librarian 

ESL,  Inc. 

ATTN 

T.  Potemra 

ATTN 

J.  Marshall 

ATTN 

B.  Wise 

ATTN 

C.  Prettie 

ATTN 

T.  Evans 

ATTN:  J.  Roberts 

ATTN 

J.  Newland 

ATTN 

P.  Komiske 

Ford  Aerospace  &  Communications  Corp. 

ATTN 

J.  Mattingley 

Kaman  Sciences  Corp. 

ATTN 

T.  Meagher 

General  Electric  Co. 

Space  Division 

Linkabit  Corp. 

Valley  Forge  Space  Center 

ATTN 

I.  Jacobs 

ATTN 

M.  Bortner 

Litton  Systems,  Inc. 

General  Electric  Co. 

Amecom  Division 

Re-Entry  & 

Environmental  Systems  Div. 

ATTN 

R.  Grasty 

ATTN 

A.  Steinmayer 

ATTN 

C.  Zierdt 

Lockheed  Missiles  &  Space  Co.,  Inc. 

ATTN 

S.  Lipson 

ATTN 

D.  Churchill 

ATTN 

Dept.  60-12 

General  Electric  Co. 

ATTN 

F.  Reibert 

Lockheed  Missiles  and  Space  Co.,  Inc. 

ATTN 

R.  Johnson 

General  Electric  Co. -TEMPO 

ATTN 

W.  Imhof 

ATTN 

T.  Stevens 

ATTN 

M.  Walt 

ATTN 

M.  Stanton 

ATTN 

DASIAC 

M.I.T.  Lincoln  Lab. 

ATTN 

W.  Knapp 

ATTN 

D.  Towle 

ATTN 

D.  Chandler 

ATTN 

L.  Loughlin 

General  Electric  Tech.  Services  Co.,  Inc. 

McDonnell  Douglas  Corp. 

ATTN 

G.  Millman 

ATTN 

G.  Mroz 

ATTN 

W.  Olson 

General  Research  Corp. 

ATTN 

N.  Harris 

ATTN 

J.  Garbarino 

ATTN 

J.  Moule 

ATTN 

J.  Ise,  Jr. 

Mission  Research  Corp. 

GTE  Syl vania,  Inc. 

ATTN 

R.  Bogusch 

Electronics  Systems  Grp-Eastern  Div. 

ATTN 

S.  Gutsche 

ATTN 

M.  Cross 

ATTN 

D.  Sowle 

ATTN 

R.  Hendrick 

HSS,  Inc. 

ATTN 

M.  Scheibe 

ATTN 

D.  Hansen 

ATTN 

F.  Fajen 

IBM  Corp. 

Mitre  Corp 

Federal  Systems  Division 

ATTN 

G.  Harding 

ATTN 

F.  Ricci 

ATTN 

C.  Callahan 

ATTN 

A.  Kymmel 

University  of  Illinois 

ATTN 

K.  Yeh 

Mitre  Corp. 

ATTN 

M.  Horrocks 

Institute  for  Defense  Analyses 

ATTN 

W.  Foster 

ATTN 

J.  Bengston 

ATTN 

W.  Hall 

ATTN 

E.  Bauer 

ATTN 

J.  Aein 

Pacific-Sierra  Research  Corp. 

ATTN 

H.  Wolfhard 

ATTN 

E.  Field,  Jr. 

International  Tel.  &  Telegraph  Corp. 

Pennsylvania  State  University 

ATTN 

G.  Wetmore 

Ionosphere  Research  Lab. 

ATTN 

Technical  Library 

ATTN 

Ionospheric  Research  Lab. 

Jaycor 

Photometries,  Inc. 

ATTN:  S.  Goldman 

ATTN:  I.  Kofsky 

Jaycor 

Physical  Dynamics,  Inc. 

ATTN 

D.  Carlos 

ATTN 

E.  Fremouw 
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R. 

Lelevier 

ATTN:  C.  Smith 

M. 

Gantsweg 

ATTN:  D.  Hamlin 

H. 

Ory 

ATTN:  D.  Sachs 

R. 

Turco 

ATTN:  J.  McDougall 

C. 

Griefinger 

ATTN:  E.  Straker 

w. 

Wright,  Jr. 

ATTN:  L.  Linson 
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R  &  D  Associates 

ATTN:  R.  Lelevier 
ATTN:  M.  Gantsweg 
ATTN:  H.  Ory 
ATTN:  R.  Torco 
ATTN:  C.  Griefinger 
ATTN:  U.  Wright,  Jr. 

ATTN:  B.  Gabbard 
ATTN:  W.  Karzas 
ATTN:  C.  MacDonald 
ATTN:  F.  Gilmore 

R  &  D  Associates 

ATTN:  L.  Delaney 

Rand  Corp. 

ATTN:  C.  Crain 
ATTN:  E.  Bedrozian 

Riverside  Research  Institute 
ATTN:  V.  Trapani 

Rockwell  International  Corp, 

Collins  Telecommunications  Sys.  Division 
ATTN:  J.  Kristof 

Santa  Fe  Corp. 

ATTN:  E.  Ortlieb 

Utah  State  University 
Space  Measurements  Lab. 

ATTN:  L.  Jensen 
ATTN:  K.  Baker 

Visidyne,  Inc. 

ATTN:  J.  Carpenter 
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Science  Applications,  Inc. 
ATTN:  D.  Divis 

Science  Applications,  Inc. 
ATTN:  SZ 

SRI  International 

ATTN:  A.  Burns 
ATTN:  G.  Price 
ATTN:  R.  Leadabrand 
ATTN:  M.  Baron 
ATTN:  G.  Carpenter 
ATTN:  R.  Livingston 
ATTN:  C.  Rino 
ATTN:  W.  Jaye 
ATTN:  W.  Chesnut 
ATTN:  D.  Neil  son 
ATTN:  G.  Smith 

Teledyne  Brown  Engineering 
ATTN:  R.  Deliberis 

Tri-Corn,  Inc. 

ATTN:  D.  Murray 

TRW  Defense  &  Space  Sys.  Group 
ATTN:  S.  Altschuler 
ATTN:  R.  Plebuch 
ATTN:  D.  Dee 


